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We report muon spin rotation (/iSR) measurements of single crystal Ba(Fei_i:Co2,)2As2 and 
Sr(Fei_a;Coa;)2As2. From measurements of the magnetic field penetration depth A we find that 
for optimally- and over-doped samples, 1/A(r — ^ 0)^ varies monotonically with the superconducting 
transition temperature Tc- Within the superconducting state we observe a positive shift in the 
muon precession signal, likely indicating that the applied field induces an internal magnetic field. 
The size of the induced field decreases with increasing doping but is present for all Co concentrations 
studied. 

PACS numbers: 76.75.+i 



I. INTRODUCTION 

Of the various families of iron pnictide supercon- 
ductors, the so-called 122 family has been extensively 
studied due to their high Tc's and the abihty to grow 
single crystals. This family includes BaFe2As2 and 
SrFe2As2. Unlike the cuprates, these materials are 
quite robust against in-plane disorder, brought about 
by electron-doping for Fe atoms either by Co, Ni or 
other transition metals. The transition temperatures 
remain fairly high for these substitutions, with Tc 
= 22K for Ba(Feo.926Coo.o74)2As2lTl |g, 20.5K for 
Ba(Feo.952Nio.o48)2As2l3], 23K for Ba(Feo.9Pto.i)2As2|l, 
14K for Ba(Feo.96iRho.o39)2As2[5], 19.5K 

for Sr(Feo.8Coo.2)2As2[6] and 9.5K for 

Sr(Feo.925Nio.o75)2As2l7]. 

Measurements of the penetration depth and superfluid 
density have attempted to address the nature of the su- 
perconducting gap symmetry. NMR has shown the lack 
of a coherence peakfSJ fS|, indicative of unconventional 
pairing. Similarly, tunnel-diode resonator measurements 
in Ba(Fei_2,Co3;)2As2 also show power-law temperature 
dependences for the penetration depthfTTHfTTl, which are 
interpreted in terms of gap nodes. Other measurements 
see a constant superfluid density at low temperatures [12 
[13], indicating an s-wave gap. Likewise, the possibility 
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of multi-band superconductivity in the pnictides has also 
been studied by analyzing the temperature-dependence 
of the superfluid density [T4l[T5]. 



II. EXPERIMENTAL 

Muon spin rotation (/iSR) is a powerful local micro- 
scopic tool for characterizing the magnetic properties of 
materials, in superconducting or other states. A thor- 
ough description of the application of /iSR to studies of 
superconductivity can be found elsewhere ^6j. In a trans- 
verse field (TF) /xSR experiment, spin-polarized positive 
muons are implanted one at a time into a sample. Each 
muon spin precesses around the local magnetic field until 
the muon decays into a positron, which is preferentially 
ejected along the direction of the muon spin at the time of 
decay (as well as two neutrinos which are not detected). 
In the presence of a vortex lattice, the spatial variation 
of the magnetic field distribution results in a dephasing 
of the muon spin polarization and a relaxation of the 
precession signal. A Fourier transform of the spin polar- 
ization function essentially reveals the field distribution 
which exhibits a characteristic Abrikosov lineshape. The 
lineshape (or equivalently the relaxation function in the 
time domain) depends on the lattice geometry, magnetic 
field penetration depth A, coherence length, ^, and the 
amount of lattice disorder. As a result, careful analy- 
sis of the relaxation function allows these microscopic 
parameters to be determined in the vortex state. Such 
measurements demonstrated the presence of gap nodes 



characteristic of d-wave pairing in high quahty single 
crystals of YBa2Cu306.97[16]. In ceramic samples this 
anisotropic lineshape is generally not observed, rather 
the broadened line is generally well described by a Gaus- 
sian distribution; however, the width of this distribution 
(the Gaussian relaxation rate) a has been shown to be 
proportional to the superfluid density divided by the ef- 
fective mass a ex ns/m* ex 1/A^|T71[TH]. Previous studies 
of cuprates found that extrinsic effects in ceramics can 
result in the correct temperature dependence of the su- 
perfluid density being masked; for this reason, reliable 
measurements of the superfluid density require the use 
of single crystals and the observation of an anisotropic 
lineshape characteristic of a vortex lattice. 

High quality single crystals of Ba(Fei_2,Coa;)2As2 with 
x=0.061, 0.074, 0.107 and 0.114 were grown at Ames 
from self flux as described in detail elsewhere [2] • Some 
measurements of the x = 0.074 sample were reported 
previously 1 14i. A single crystal of Sr(Fei_2;Co2,)2As2 
with X = 0.13 was grown at Maryland, also from self 
flux'19]. The crystals, each of roughly 1 cm^ area, were 
mounted in a helium gas flow cryostat on the M20 sur- 
face muon beamline at TRIUMF, using a low background 
arrangement such that only positrons originating from 
muons landing in the specimens were collected in the ex- 
perimental spectra. Zero field /iSR measurements of each 
sample confirmed that no magnetic order or spin freezing 
was present in any of the samples. 



III. PENETRATION DEPTH 



Fourier transforms of the transverse field (TF)-/iSR 
spectra (a representative set are shown in Fig. n]) exhibit 
the anisotropic lineshape characteristic of an Abrikosov 
lattice, indicating the presence of a least locally well- 
ordered vortices in the superconducting state. All of the 
the Fourier transform lineshapes are consistent with a tri- 
angular vortex lattice; for example, a square lattice would 
have the frequency corresponding to the most likely field 
(the peak of the lineshape) much more separated from 
the minimum field in the field distribution. We analyzed 
the data by fitting the spectra to an analytical Ginzburg- 
Landau model which allows us to calculate theoretical 
/iSR time spectra as a function of the vortex lattice ge- 
ometry, magnetic field penetration depth (A) and coher- 
ence length (f). We included the effects of vortex lat- 
tice disorder in our analysis via an additional Gaussian 
broadening of our /iSR spectrum[17, 20 , where we as- 
sumed that this broadening was proportional to 1/A^ as 
observed in previous studies of cuprates and other high 
K superconductors |16j . The errors quoted in various fit 
parameters included the correlations between the vari- 
ous parameters. The fit parameters were fairly weakly 
correlated since the effect of each parameter on the re- 
laxation function is reasonably unique: the penetration 
depth affects the overall linewidth, the coherence length 
affects the high field cutoff while disorder gives an overall 



broadening of the various van Hove singularities in the 
lineshape. Consistent with our previous measurements 
of Ba(Feo.926Coo. 074)2 As2 [H], we found that the rms 
deviation of the vortex positions {{s'^)^^'^) relative to the 
vortex separation was greatest in lower fields (up to 30% 
in 0.02 T at low temperature) and smallest at the highest 
fields (about 2% in 0.1 T) and decreased with increasing 
temperature. The disorder was greatest for the samples 
with the highest Tc- 
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FIG. 1: (Color online) Fourier transforms of TF-/1SR spectra 
for Ba(Fei_a;Coi:)2As2, showing anisotropic lineshapes char- 
acteristic of an Abrikosov vortex lattice. 

Results of this analysis for 1/A^ are shown in Fig.p^for 
applied fields of 0.1 T, and 0.02 T. In conventional weak- 
coupling BCS theory, the low-temperature behaviour of 
1/A^ should be exponentially flat while the presence of 
gap nodes would be reflected in low-temperature power 
law behaviour. We see in Fig.[2]that the low-temperature 
behaviour varies more rapidly than standard BCS pre- 
dictions, and also note that recent specific heat measure- 
ments have observed the possibility of gap nodes [21] . Fol- 
lowing our earlier work on the Ba(Feo. 926C00. 074)2 As2[T3] 
we fit the superfiuid density to a phenomenological two- 



gap model[221 123] which has been employed in previ- 
ous /xSR studies of LaFeAs(0,F), Ca(Fe,Co)AsO, and 
(Ba,K)Fe2As2l21], 

ns{T)^nsiO)-w6nsiAi,T)-il-w)-Sns{A2,T) (1) 

where w is the relative weight for the first gap, Ai. Here, 
the gap functions are given by. 



dn{A,T) 



2n,(0) 
ksT JO 



f{e,T)-[l-.f{e,T)]de (2) 



where f{e,T) is the Fermi distribution given by, 

/(e, T) = (1 + eV^-^+^mV^-«^)-i (3) 

Here, A^ (i = 1 and 2) arc the energy gaps at T = 0, and 
Ai(T) were taken to follow the standard weak-coupled 
BCS temperature dependence. This model reduces to a 
single-gap BCS model when w = 1. The size of the gaps, 
Ai and A2, and Tq were fit globally, while ris(O) and the 
weighting factor, w, were allowed to be field-dependent. 
The fit values of 1/A^ are shown by the fit lines on Fig.p] 
We see that this two gap model fits the observed temper- 
ature dependence of 1/A^ for each of the samples and 
fields measured. Single gap fits did not give satisfac- 
tory results (when the BCS gap function was used). For 
most of the samples we obtained the larger gap value 
2A//cbTc ~ 3.7 which is close to the weak coupled BCS 
value of 3.5. For the Sr(Feo.87Coo.i3)2As2 the larger gap 
was 2A/kBTc « 2.7, less than the BCS value. For the 
Ba(Fei_a;Coa;)2As2 samples with x = 0.107, 0.114 most 
of the weight was on the smaller gap. The results for 
these three samples give a stronger low temperature de- 
pendence to the superfluid density than for a single-gap 
weak-coupled BCS system, as can be seen on Fig. [2] This 
steeper temperature dependence may possibly reflect a 
non-s-wave gap as has been interpreted by tunnel diode 
resonator measurements. [TUl [TT] We don't have enough 
data points, especially at temperatures below 2 K, to 
make a definitive statement regarding the presence of 
gap nodes. We are able, however, to reliably extrapolate 
the superfluid density to obtain a good estimate of the 
magnitude of 1/A^(T — )• 0). Examining the behaviour of 
1/A^ for the different samples in Fig. p^ we see that there 
is considerable variation in 1/A^(T — >■ 0). Over the range 
of dopings and fields studied, the value of 1/A^(T -> 0) 
varies from 5 ^m~^ to nearly 30 /irn"^, more than half 
an order of magnitude. 

There is considerable field dependence in 1/A^(T -^ 0) 
for the Ba(Fei_a;Coa:)2As2 samples with x = 0.061 and 
X = 0.074 which is essentially absent for the higher 
doped samples with the smaller superfluid density. We 
first noted this large field dependence in our study of 
Ba(Feo.926Coo.o74)2As2[T3]; subsequent studies on other 
pnictides have seen similar behaviour. |25) The field de- 
pendence in the density of states in a multiband super- 
conductor has been calculated by Ichioka et a/.[2S] who 
noted that a strong field dependence is expected for fields 



Tc 



Ba(Feo.939Coo.o6i)2As2 23.6 

Ba(Feo.926Coo.o74)2As2 22.1 

Ba(Feo.899Coo.ioi)2As2 14.1 

Ba(Feo.89Coo.ii)2As2 10.3 

Sr(Feo.87Coo.i3)2As2 16.2 



Ao(0.02r) 



189.4 ±1.1 
224.2 ±0.6 
332.2 ±2.2 
453.8 ±2.6 

325.5 ±0.5 



Ao(O.ir) Ao(B = 0) 



240.5 ±2.0 182.6 ±1.4 
277.4 ±1.0 216.8 ±0.7 
348.3 ±4.6 329.3 ±3.4 
448.0 ±2.4 454.9 ±3.6 
339.8 ±0.6 322.8 ±0.7 



TABLE I: Results of fitting 1/A^(r) to Eqn.[l]for Tc(K), Ao 
(nm) in fields of 0.02 T and 0.1 T. Also shown are values of 
Ao(nni) extrapolated to zero field. 



on the order of the smaller gap size. Results of previous 
/iSR measurements of a variety of multi-gap supercon- 
ductors are described in Ref. [H]. In those materials 
(such as NbSe2) the origin of the field dependence is the 
loosely bound core states associated with the smaller gap. 
With increasing field these states become more delocal- 
ized and affect the field distribution seen by the muon 
ensemble. The field dependence we observe in this study 
is larger than we would expect to be due to the smaller 
gap A2 and may possibly indicate that an anisotropic gap 
(perhaps with nodes) might be more appropriate for the 
smaller gap than the uniform gap model used to fit the 
temperature dependence of the superfluid density. We 
note that different gap symmetries on different parts of 
the Fermi surface might resolve the apparent discrepan- 
cies between different techniques that probe the normal 
state carrier concentration (eg. tunnel diode oscillator, 
microwave) and superfluid carrier concentrations (/iSR). 
To estimate the zero applied field values of Aq we have 
performed a linear extrapolation of the fit values of 1/Aq 
measured in 0.02 T and 0.1 T and included the resulting 
Ao(S -^ 0) values in Table |!) 
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FIG. 2: (Color online) Measured 1/A^ for Ba(Fei_^Co^)2As2 
and Sr(Fei_^Co^)2As2 measured in TF=0.02 T (filled sym- 
bols) and 0.1 T (open symbols). 

Fig. pi shows the extrapolated values of 1/A^(T — > 0) 
and the fit values of Tc as a function of the level of Co 
doping X for Ba(Fei_a;Co2;)2As2 and Sr(Fei_a;Co2;)2As2. 



We see that above x = 0.06 in Ba(Fei_a;Coa:)2As2 the su- 
perconducting Tc decreases with increasing Co substitu- 
tion, in agreement with previous work. Additionally, the 
Tc for Sr(Fei_a;Co2:)2As2 does not lie on the same curve 
as for the Ba(Fei_2:Coa;)2As2 family; the location of the 
superconducting phase dome within the phase diagram 
is different for the two families. The lower panel of Fig. [3] 
shows the evolution of the extrapolated 1/A^(T -> 0). 
Within the Ba(Fei_2;Coa;)2As2 family, there is a mono- 
tonic decrease in 1/A^(T —>■ 0) and again, the point for 
Sr(Fei_a;COa;)2As2 does not he on the same curve. 
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FIG. 3: (Color online) Superconducting Tc's and 1/A^(T -i> 
0) for Ba(Fei_a;Coi)2As2 and Sr(Fei_a;Coi)2As2 as a function 
of Co concentration x, measured in TF=0.02T and O.IT, and 
extrapolated to B = 0. The open points and dashed lines are 
the measured Tc's and the superconducting dome taken from 
[2] for Ba(Fe,Co)2As2 and [3 [17] for Sr(Fe,Co)2As2. 

Muon spin rotation measurements on a wide variety 
of cuprate and other exotic superconductors have re- 
vealed a strong, roughly linear correlation between the 
superconducting transition temperature and the extrap- 
olated zero temperature superfluid density divided by 
the effective mass. [28 This relation is not expected in 
standard BCS theory, implying that a different mecha- 
nism is responsible for superconductivity in these sys- 
tems. We plot our fit values of Tc vs. 1/A2(T -> 0) 
in Fig[4]for Ba(Fei_3,Co3;)2As2 and Sr(Feo.87Con.i3As2)2. 
In contrast to the plots on Fig. [3] the points for all 
of the samples lie close to common curves for both 
20 niT and 100 mT. We see that in single crystals of 
Ba(Fei„a,C0j;)2As2 and Sr(Fei_2.Co2,)2As2, the super- 
conducting Tc is apparently determined by the carrier 
density divided by the effective mass. 

Specific heat measurements [55] of the superconduct- 
ing transition found that the specific heat jump at Tc 
divided by Tc was correlated with Tc as ACp/Tc oc T^. 
Our results for 1/A(T — >■ 0)^ are shown, along with spe- 
cific heat jump for Sr(Feo.87Coo.i3)2As2 and the results 
of Bud'ko et oL[3S] in Fig. [s] In agreement with the 
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FIG. 4: (Color online) Superconducting Tc vs. 1/A^(T -> 0) 
for Ba(Fei_a;Coi,)2As2 and Sr(Fei_i,Coa;)2As2 as a function 
of Co concentration x, measured in TF=0.02T and O.IT and 
extrapolated to B = 0. 



specific heat, we find that 1/A(T —> 0)^ can be well de- 
scribed by a straight line with slope n « 2 as indicated 
by the dashed line. The common variation of the su- 
perfiuid density [36] and the specific heat jump ACp/Tc 
and Tc,[37j as a function of carrier doping was first noted 
in overdoped Tl2Ba2CuOg+5 cuprates[38". The present 
case of the FeAs superconductors, shown in Fig. [5J ex- 
hibits commonalities to the cuprates in this regard. 
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FIG. 5; (Color online) 1/A^(T -> 0) vs. superconducting Tc 
for Ba(Fei_^Coi,)2As2 and Sr(Fei_i,Coa;)2As2 as a function 
of Co concentration a;, measured in TF=0.02T and O.IT and 
extrapolated to B = 0, plotted with ACp/Tc from Ref. [35] . 
The dashed line has a slope n = 2 and is a guide to the eye. 

There are two contributions to the measured pene- 
tration depth with chemical substitution. First, doping 
changes the carrier concentration which directly changes 



the London penetration depth Al via I/Xj^ ex ns/m*. 
If a system is in the dirty Umit, the measured pene- 
tration depth is actually an effective penetration depth 
^eff — ^l{^ + Co/0' where ^o is the coherence length 
and I is the mean free path. [2^ The upper critical field 
is quite large in these systems, taking Hc2 « 50 T gives 
an estimate of ^o = 2.5 nm. In order to estimate the 
mean free path, a reasonable value of the Fermi velocity 
is needed. Due to the unclear situation of the nature of 
the Fermi surface in the pnictides, an estimate of the pair- 
breaking effect is unlikely to be accurate. Optical conduc- 
tivity measurements have directly detected the opening 
of the superconducting gap in Ba(Fei_^Coa;)2As2 with 
X = 0.1[3n], X = 0.07[3l] and x = 0.065[32 . Broadening 
of the normal-state zero frequency Drude conductivity 
indicates the presence of significant normal-state scatter- 
ing. Some authors have argued that strong pair-breaking 
can account for a number of effects, including the specific 
heat jump at Tc and the behaviour of dHc2/dT[S5 . As- 
suming strong pair-breaking, Kogan[33] has found that 
1/A(0) ex Tc, in agreement with our results in Fig. ^ 
However, the existence of such strong pair breaking is 
not yet proven. Although scattering is clearly present 
in these systems, it is unlikely that reasonably modest 
changes in the dopant concentration (of a few percent) 
could cause such a dramatic change in the scattering so 
as to dominate the penetration depth and as such, sub- 
stantial changes in the superfluid density are apparently 
occurring with chemical substitution. 

Hall effect measurements show that the normal-state 
carrier concentration increases monotonically with in- 
creased chemical substitution 39]. If in fact the super- 
fluid density decreases with increasing doping above the 
maximum Tc, then this implies that not all of the carri- 
ers join the condensate below Tc. This segregation into 
superconducting and normal fluids could be in reciprocal 
space, if superconductivity occurs on only some parts of 
the Fermi surface. This could also occur in real space 
with phase separation into normal and superconduct- 
ing regions. Previous /iSR measurements of overdoped 
Tl2Ba2CuO6+5[3nill0] exhibited similar behaviour, with 
increased normal state doping and a loss of superconduct- 
ing carrier density. Real space phase separation has been 
seen in other /iSR measurements of both hole- |3TJ 35] and 
electron-doped 03] pnictide superconductors. 

Phase separation (either in real space or reciprocal 
space) should leave a residual normal fluid whose spectral 
weight should increase with Co substitution and which 
should be apparent in measurements of optical conduc- 
tivity. Recent optical measurements by Gorshunov et 
al, in Ba(Fei_j;Coa;)2As2 with x — 0.1 do in fact show 
appreciable residual conductivity well within the super- 
conducting state which may be evidence of this resid- 
ual normal fluid |30|. Additional measurements at low 
frequencies for a range of doping are needed to further 
test this hypothesis. Gofryk et al. reported specific heat 
measurements of Ba(Fei_a;COa;)2As2 for a range of Co 
concentrations. [44 . They found a substantial normal 



fluid response (residual linear specific heat contribution) 
which increased with Co concentration for x > 0.08. 

Our /zSR results in Co substituted BaFe2As2 and 
SrFe2As2 indicate that the vortex lattice exists through- 
out the samples which indicates that any phase separa- 
tion is cither in real space with a characteristic length- 
scale much less than the penetration depth (perhaps the 
coherence length) or in reciprocal space. A model for 
real space phase separation for overdoped cuprates has 
been discussed in Ref. |3S] . If the phase separation oc- 
curs in momentum space, it could originate perhaps from 
only some of the multiple bands in these systems par- 
ticipating in the pairing. Angle-resolved photoemission 
measurements ^5] have shown that above x ~ 0.08 over- 
doped electrons fill the hole Fermi surface at the Brillouin 
zone centre, resulting in a loss of interband scattering. If 
this scattering is responsible for pairing, then the loss 
of the hole states with substitution could reduce the su- 
perfluid density, even though the normal state carrier 
concentration increases with doping. 



IV. PARAMAGNETIC FREQUENCY SHIFT 

When fitting the /nSR time spectra to our analyti- 
cal Ginzburg-Landau model, one of the fitted parame- 
ters is the average muon precession frequency v^,. In 
the normal state, this precession frequency is given by 
i^p, = {I + Kf^)j^Bc:>^t where Boxt is the externally ap- 
plied magnetic field, 7^ is the muon gyromagnetic ra- 
tio and K^ is the Knight shift. In the superconducting 
state the muon precession frequency is generally slightly 
reduced from the normal state value due to fiux expul- 
sion; for thin plate-like samples this reduction is gener- 
ally quite small due to the demagnetizing factor. We 
show the fitted values of the fractional shift in the pre- 
cession frequency relative to its normal state value i'^tv in 
Fig. [6J We see that except for a negative shift right below 
T(7 for some samples due to bulk screening, all samples 
have increasing frequency shifts with decreasing tempera- 
ture in the superconducting state. A similar positive fre- 
quency shift has also been reported by Khasanov et al. in 
SrFei.75Coo.25As2[17]. A positive value of f^/f^Ar — 1 in- 
dicates that the field at the muon site is actually greater 
than the applied field. Since bulk screening can only 
contribute a negative frequency shift, we need to find a 
different explanation for our observed positive shifts. 

The fractional shift within the superconducting state 
is considerably larger in the 0.02T data than in the O.IT 
runs as shown in Fig.[6j In fact, the absolute value of the 
shifts (i^p — I'fj.N) is roughly the same for the two fields. 
The shifts are also largest for the samples with the high- 
est Tc and highest superfluid density ris/m* ex 1/A^. 
Previous fi'SR studies of the electron-doped cuprate su- 
perconductor Pr2_3;Cea;Cu04 also exhibited a positive 
frequency shift below Tc which was interpreted as ev- 
idence of field-induced magnetism 48 . In that case, the 
absolute shift decreased with increasing field (not just 
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FIG. 6: (Color online) Fractional shift of muon precession 
frequency v^/v^n — 1 relative to the normal state frequency 
fpAT. Superconducting transition temperatures indicated by 
triangle symbols. 



the fractional shift), indicating that the induced fields 
were perpendicular to the applied field. In the present 
case, the fact that the absolute shift is roughly field- 
independent indicates that the induced moments must 
be parallel to the applied field and have a ferromagnetic 
character (antiferromagnetic fields would split the pre- 
cession line, rather than shift it). We note that these 
field-induced ferromagnetic fields would not be apparent 
in bulk susceptibility measurements, since they would be 
screened by supercurrents on the surface of the sample. 
In each sample the paramagnetic frequency shift sets in 
at the superconducting Tj, of each particular sample, im- 
plying that it is a property of the superconducting state. 
One possible source of such a field could be a spin triplet 
pair state where the Cooper pairs possess a non-zero an- 
gular momentum. However, other explanations are also 
possible and further experiments will be required to de- 



termine the microscopic origin of these internal fields. 



V. CONCLUSIONS 

We have measured the London penetration depth in 
single crystals of Ba(Fei„a;Co2;)2As2 (withx > 0.061)and 
Sr(Feo.87Coo.i3)2As2 using muon spin rotation. The tem- 
perature dependence of 1/A^ can be fit by a two-band 
model where the gaps follow the BCS temperature de- 
pendence. In the more highly-doped samples we find that 
dominant gap magnitude is considerably smaller than the 
weak-coupled BCS result, implying that the gap must be 
highly anisotropic and could possess nodes in this dop- 
ing regime. Our results demonstrate that I/A^(T — > 0) 
varies roughly quadradically with the superconducting 
transition temperature Tc. We find that the superfluid 
density divided by the effective mass n^/m* ex I/A* de- 
creases as normal state charge carriers are added. This 
implies that a form of electronic phase separation (either 
in real or reciprocal space) occurs in these systems. We 
observe a paramagnetic frequency shift in all specimens 
below Tc, the magnitude of which is roughly independent 
of field but decreases with increasing doping. 
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